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Abstract — In this paper we present several methods and con- 
structions to generate codes for correction of a multidimensional 
cluster-error. The goal is to correct a cluster-error whose shape 
can be a box-error, a Lee sphere error, or an error with an 
arbitrary shape. Our codes have very low redundancy, close to 
optimal, and large range of parameters of arrays and clusters. 
Our main results are summarized as follows: 

1) A construction of two-dimensional codes capable to cor- 
rect a rectangular-error with considerably more flexible 
parameters from previously known constructions. Another 
advantage of this construction over previously known con- 
structions is that it is easily generalized for D dimensions. 

2) A novel method based on D colorings of the D-dimensional 
space for constructing D-dimensional codes correcting D- 
dimensional cluster-error of various shapes. This method is 
applied efficiently to correct a D-dimensional cluster error 
with parameters not covered efficiently by the previous 
constructions. 

3) A transformation of the D-dimensional space into another 
-D-dimensional space in a way that a D-dimensional Lee 
sphere is transformed into a shape located in a D- 
dimensional box of a relatively small size. This transforma- 
tion enables us to use the previous constructions to correct 
a D-dimensional error whose shape is a D-dimensional Lee 
sphere. 

4) Applying the coloring method to correct more efficiently a 
two-dimensional error whose shape is a Lee sphere. The 
D-dimensional case is also discussed. 

5) A construction of one-dimensional codes capable to correct 
a burst-error of length b in which the number of erroneous 
positions is relatively small compared to b. This construc- 
tion is generalized for D-dimensional codes. 

6) Applying the construction for correction of a Lee sphere 
error and the construction for correction of a cluster-error 
with small number of erroneous positions, to correct a D- 
dimensional arbitrary cluster-error. 

Index Terms — burst-error, burst-locator code, cluster- 
correcting code, coloring, Lee sphere, multidimensional code. 



I. Introduction 

In current memory devices for advanced storage systems 
the information is stored in two or more dimensions. In 
such systems errors usually take the form of multidimensional 
bursts. Usually, a cluster of errors either will be affected by 
the position in which the error event occurred or will be of an 
arbitrary shape. But, since an arbitrary cluster-error is hard to 
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correct efficiently it is common to assume some type of cluster- 
error (as any arbitrary cluster is located inside a cluster with 
a certain shape, e.g., any two-dimensional cluster is located 
inside a rectangle). These types of errors can be of specific 
shapes like rectangles or Lee spheres. We will consider these 
types of errors as well as arbitrary cluster-errors. The main 
measure to compute the efficiency of a cluster-error correcting 
code is its redundancy. If we want to design a code which 
corrects one multidimensional cluster-error with volume B 
(of an arbitrary or a specific shape) then the redundancy of 
the code r satisfies r > 2B. This bound, known as the 
Reiger bound [14], is attained for binary two-dimensional 
codes, which correct a rectangular error, constructed recently 
by Boyarinov [7]. If the volume of the array is N then the 
redundancy of the code must also satisfy r > \og 2 N + B — 1 
(usually r— |~log 2 7V] > B). The difference r— |~log 2 7V] will be 
called the excess redundancy of the code [1], [4] (even so our 
definition is slightly different). Abdel-Ghaffar [1] constructed 
a binary two-dimensional code which corrects a burst with a 
rectangle shape for which r = \log 2 N~\ + B. The code has a 
few disadvantages: very limited size, complicated construction, 
and there is no obvious generalization for higher dimensions. 

Our goal is to design codes which are capable to correct a 
cluster-error whose shape is a box, a Lee sphere, or an arbitrary 
shape. The method should be able to work on two-dimensional 
codes and multidimensional codes, and the parameters of the 
size of the codewords and the size of the cluster are as 
flexible as possible. There will be a price for our flexibility 
and our ability to generalize a two-dimensional construction 
into multidimensional construction. This price will be in the 
excess redundancy. While the two-dimensional codes which 
correct a rectangle-error of Abdel-Ghaffar [1] have optimal 
excess redundancy, the excess redundancy of our codes is 
only close to optimality. Moreover, the novel methods enable 
us to correct a cluster whose shape is a Lee sphere and an 
arbitrary cluster, with excess redundancy close to optimal or 
very low, depending on the exact parameters. Previous to our 
methods the way to correct such a cluster-error was to use a 
code which corrects a box-error in which the cluster-error is 
located, a method for which the excess redundancy is far from 
optimality. 

The rest of the paper is organized as follows. In Section HIl 
we briefly survey some of the known constructions which are 
essential to understand our results. In Section [HI] we present 
a construction for codes which correct a multidimensional 
box-error. The construction is a generalization and a mod- 
ification of the construction of Breitbach, Bossert, Zybalov, 
and Sidorenko [8] for correction of bursts of size b\ x 6 2 . 
Better codes are constructed when the volume of the D- 
dimensional box-error is an odd integer. These constructions 
and the constructions which follow use auxiliary codes, called 
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component codes, one code for each dimension. In Sec- 
tion [IV] we present a novel method for correction of a D- 
dimensional cluster. The construction uses D colorings of the 
£>-dimensional space. The construction of Section [III] is a 
special case of this construction. The new construction enables 
us to handle different burst patterns. In Section [IV] we use 
this method to handle /^-dimensional box-error, where the 
volume of the box is an even integer. In Section [V] we discuss 
how to correct a D-dimensional cluster-error whose shape is 
a Lee sphere. Two types of constructions are used. The first 
one uses a transformation of the D-dimensional space into 
another Z?-dimensional space in a way that each Lee sphere is 
transformed into a shape located inside a reasonably small D- 
dimensional box, so that we can use the constructions of the 
previous sections. The transformation is especially efficient 
for two-dimensions. The second construction uses colorings 
as in Section [IV] For two-dimensional array the colorings 
that we use result in codes with excess redundancy close to 
optimality (where our measure for optimality is the lower 
bound on the excess redundancy), which improves on the 
construction obtained by the two-dimensional transformation. 
The generalization for multidimensional Lee sphere errors 
usually does not make the same improvement. This is also 
discussed in Section[V] In Section[VT]we show how we handle 
bursts of size b, where the number of erroneous positions is 
limited. First, we present a construction for one-dimensional 
codes and afterwards we generalize it into D-dimensional 
codes. In Section [VTIl we combine the constructions of codes 
which are capable to correct Lee sphere error and the con- 
struction capable to correct a burst with a limited number 
of erroneous positions for a construction of codes capable to 
correct arbitrary bursts. In Section IVIIII we describe codes with 
the same or slightly better parameters than the parameters of 
the codes from the previous sections by using parity-check 
matrices. Finally, a conclusion and a list of problems for 
further research are given in Section [IX] 

II. Known Constructions 

Five constructions are important to understand our construc- 
tion and their comparison with previous results. 

• Abdel-Ghaffar, McEliece, Odlyzko, and van Tilborg [3] 
construction of optimum binary cyclic burst-correcting 
codes. 

> Abdel-Ghaffar [2] construction of optimum cyclic burst- 
correcting codes over GF(<7). 

• Abdel-Ghaffar construction [2] of two-dimensional codes 
which correct rectangular-error of size b\ x fe 2 with excess 
redundancy 6162. 

« Breitbach, Bossert, Zybalov, and Sidorenko construc- 
tion [8] of two-dimensional codes for correction of a 
(61 x &2)-rectangular-error by using vertical and horizontal 
component codes. 

> Abdel-Ghaffar, McEliece, and van Tilborg construc- 
tion [4] of two-dimensional burst identification codes, 
which are used to identify the shape of an error and 
together with burst location codes are used for correction 
of a two-dimensional cluster. 



Most of two-dimensional codes known in the literature are 
designed to correct a single cluster-error of size b\ x £> 2 [LJ, 
[4], [7], [8], [11] (only in some recent papers [6], [9], [17] it 
is assumed that the cluster-error can have an arbitrary shape). 
Two of these methods are important in our discussion. Abdel- 
Ghaffar [1] gave a construction of such n\ x n 2 code with 
excess redundancy 6162. One disadvantage of his method is 
that 77,2 must be considerably larger than n\ (with a possible 
exception when 62 < 2, subject to a list of restrictive 
conditions), and the existence of the code depends on series 
of restricted conditions. The main goal of his construction was 
to show that for any given integers b\ and 62 there exists a 
cyclic (61 x 62) -cluster-correcting code of some size ri\ x 122 
having optimal excess redundancy. Therefore, the size of the 
array was not a factor in his construction. His construction is 
a generalization of the optimum cyclic one-dimensional codes 
which correct a single cyclic burst of length b [2], [3]. Over 
GF(q) such code has length n, redundancy r, and it can correct 

1 — b+li 

a single cyclic burst of length b > 1, where n = q q _ 1 — . 
The existence of such codes was obtained by the following 
necessary and sufficient conditions. 

Theorem 1: If a polynomial g(x) generates an optimum b- 
burst-correcting code over GF(q), then it can be factored as 
g(x) — e(x)p(x), where e(x) and p(x) satisfy the conditions: 

1) e(x) is a square-free polynomial of degree b — 1 which 
is not divisible by x such that h e and m e are relatively 
primes to q — 1, where h e and m e are the period of e(x) 
and the degree of the splitting field of e(x), respectively. 

2) p(x) is an irreducible polynomial of degree m > b + 1 
and period q q ~i such that m and q — 1 are relatively 
primes and m = (mod m e ). 

A monic polynomial over GF(q) which satisfies condition 1) 
of Theorem Q] will be called a b-polynomial. 

Theorem 2: Let e(x) be a 6-polynomial over GF(q). Then, 
for all sufficiently large m relatively prime to q — 1 such that 
m = (mod m e ), where m e is the degree of the splitting 
field of e(x), there exists an irreducible polynomial p(x) of 
degree m such that e(x)p(x) generates an optimum fe-burst 
correcting code of length q q ~i ■ 

Remark: If the polynomial p(x) in Theorems[T]and[2]is binary 
then p{x) in the Theorems is a primitive polynomial. 

The second method is due to Breitbach, Bossert, Zybalov, 
and Sidorenko [8], who gave three constructions of binary two- 
dimensional codes of size n\ x n 2 which correct a rectangular- 
error of size 61 x & 2 - Their goal in presenting these construc- 
tions was not to obtain low excess redundancy, which is one of 
the goals in our constructions, but to present new constructions 
of codes with relatively large array size and redundancy close 
to the Reiger's bound. We will use ideas from one of the 
constructions which will be called Construction BBZS. 

A codeword {c^ } of the construction has size n\ X 12% = 
2 b2 x 2 bl with 46i&2 redundancy bits located in positions 
{(U) : 0<*<26i-l,na-6a<J<na-l}U{(i,i) : 
ni — bi < i < ri! — 1, < j < 2b 2 — 1} (see Fig. [T). These 
bits are set initially to be zeroes. Two temporary component 
codes are being used, a vertical code and an horizontal code 
(see Fig. [T). We will describe the construction of the vertical 
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code. We note an earlier construction [16] which use similar 
component codes. 

For each row i = 2b\, . . . , n\ — 1, 62 parity check bits are 
generated, pu, I = 0, 1, . . . , 62 — 1 is computed as 



Pit 



j=£,£+b 2 ,£+2b 2 ,...,j<n 2 



(1) 



The parity bits pu, I = 0, 1, . . . , 62 — 1 generate afterward 
a symbol p. = (piOiPiii ■ • ■ ,Pib 2 -x) from the extension field 
GF(2 b2 ). The symbols p., i = 2b\, . . . , n\ — 1 are considered 
as the information symbols of a Reed-Solomon (RS) code of 
length raj., dimension rii — 2&i, and minimum distance d = 
2b\ + 1. By the encoding procedure of the RS code we obtain 
2&i redundancy symbols p., 2 = 0, ... , 2&i — 1, and the 2&1&2 
upper right corner redundancy bits of the array are computed 
in a way that ([T|i will hold for i — 0,1, ... , 2b% — 1 and £ — 

0. 1. . . . , 62 — 1- The encoding procedure of the horizontal code 
is done in the same manner, where all the 4&1&2 redundancy 
bits of the array are assumed to be zeroes (it is possible to 
encode also with their new computed values, but we want to 
follow the construction exactly as in [8]). 

In the decoding procedure each row generates 62 parity 
bits according to ([TJ such that a word of length n\ over 
GF(2 b2 ) is received (the redundancy bits of the horizontal 
code are assumed to be zeroes). Assuming that the error 
occurred in the array can be confined inside a rectangle of 
size b\ x 62 ■ The generated word, of the vertical code, has 
at most b\ erroneous symbols, which can be corrected by 
the decoding procedure of the RS code. The same process 
is implemented for the horizontal code. The positions of the 
erroneous elements in the vertical codeword induce the rows 
in which errors occurred in the array. The positions of the 
erroneous elements in the horizontal codeword induce the 
columns in which errors occurred in the array. Hence, we 
locate the positions of the b\ x bi cluster-error in the array. 
The shape of the cluster, up to 62 horizontal cyclic shifts, is 
found by the vertical code. The shape of the cluster up to 
b\ vertical cyclic shifts, is found by the horizontal code. As 
we know the location of the cluster, we can use one of the 
component codes to identify the exact shape of the error and 
to correct it. 

Remark: Note that the vertical code cannot find errors inside 
the lower left corner redundancy bits. The horizontal code 
cannot find errors inside the upper right corner redundancy 
bits. But, these facts does not affect the decoding procedure, 

1. e., the vertical code is able to know the erroneous rows even 
if some the erroneous bits are the lower left redundancy bits. 

There is no obvious generalization to the construction of 
Abdel-Ghaffar [1] for multidimensional codes, while immedi- 
ate generalizations of Construction BBZS cannot support good 
redundancy or excess redundancy. One simple way to gener- 
alize this construction is to use the optimum burst-correcting 
codes of [2], [3] instead of the RS codes. The vertical compo- 

nent code over GF(2° 2 ) has length - — 2 b 2-i . redundancy 

r%, and it can correct a burst of length 61. The horizontal 
component code over GF(2 x ) has length ^ — ; 9 b 1 _ 1 > re- 
dundancy r 2 , and it can correct a burst of length b 2 . Instead of 
4&1&2 redundancy bits we will use ri&2 redundancy bits for the 
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Fig. 1. Construction BBZS 



vertical code and 6i7"2 redundancy bits for the horizontal code. 
The excess redundancy of this construction is 2&162 — 1 and 
the excess redundancy of its generalization for D dimensions 
is DB — 1, where B is the volume of the D-dimensional box 
error. 

Further improvements of this construction are presented 
in the next section. Henceforth we assume that if a D- 
dimensional code is discussed then D is a constant. Further- 
more, we assume b L > 1 for 1 < i < D; this assumption can 
be made since if for some j, bj = 1 then the cluster can be 
corrected as a (D — l)-dimensional cluster in a D-dimensional 
array. 



III. Construction for Multidimensional Arrays 

In this section we present our first idea for construction of 
multidimensional codes capable to correct a box-error of size 
61 x 62 X • • • x bo- First, we will present the two-dimensional 
version of the construction. We combine Construction BBZS 
with the constructions of Abdel-Ghaffar et al. [3] and Abdel- 
Ghaffar [2] to obtain codes with variety of parameters. The 
redundancy of the construction is kept relatively small as our 
horizontal code will find only the location of the error and not 
its shape. This idea is the first key of all our constructions. 
The second idea to reduce the redundancy is to use a binary 
horizontal code instead of a code over GF(2 fcl )- Finally, the 
structure of the construction makes it possible to generalize it 
to any dimension. The generalization is relatively quite simple, 
with low redundancy, and can be applied on a large range 
of parameters. One disadvantage is that the construction is 
defined for a box-error whose volume is an odd integer. To 
apply the construction on a box-error whose volume is an even 
integer, we have to increase the box-error artificially such that 
its volume will be an odd integer and the real box error will 
be located inside the artificial box error. This will cost us extra 
unnecessary redundancy. In the next section we will solve this 
problem by giving a novel construction for correction of a 
box-error whose volume is an even integer. 
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A. Two-dimensional codes 

The vertical component code of Construction BBZS finds 
the rows in which the burst occurred and the shape of the 
cluster up to a cyclic permutation of the columns. Hence, the 
work done by the horizontal code to find the shape of the 
cluster is redundant. Therefore, we want to find an horizontal 
component code that will determine only the first column of 
the cluster. More explicitly, a burst e = (eo, e\, . . . , e(, 2 _i), 
where e t G GF(2 bl ), for < i < b 2 - 1 found by the 
vertical code, can start at any column < i < n 2 — b 2 (See 
Fig. 0. However, if the first column of the cluster is i, then 
the cluster occurred is e' = (ei ,ei +i, . . . ,ei +/, 2 _i) where 
iq = i(mod b 2 ), and indices are taken modulo b 2 . 

Our new construction, in which the horizontal code only 
locates the first column of the cluster, is based on two lemmas. 
The first one is proved here only for the binary case. 

Lemma 1: If e%{x) = 1 + x + x 2 + ■ ■ ■ + x b ~ x and 6 is an 
odd integer then e 2 (x) is a 6-polynomial over GF(2). 

Proof: Clearly, e 2 (x) is not divisible by x. The derivative 
of x b — 1 over GF(2) is x h ~ x , and since g.c.d.(a; b — 1, x fc_1 ) = 
1, it follows that x b — 1 is a square-free polynomial and hence 
e 2 (x) is also square-free. Therefore, by Theorem Q] 62(21) is 
a 6-polynomial over GF(2). ■ 

There is an alternative more general version of Lemma Q] 

Lemma 2: Let 62(2:) be the polynomial e 2 (x) = 1 + x + 
x 2 + ■ ■ ■ + x fc2_1 over GF(2 bl ), where b\ and b 2 are positive 
integers. Assume that the following conditions hold: 

1) gcd(6 2 ,2) = l, 

2) gcd(6 2 ,2 fc i -1) = 1, 

3) gcd(0(6 2 ),2 bl -l) = l. 

Then, e 2 (x) is a &2-p°ly nom i a l over GF(2 bl ). 

We omit the proof of the lemma (see [18]) as the con- 
struction which uses the lemma, and is described in this 
subsection, has inferior redundancy than the one described in 
the next subsection and uses Lemma[T] We will compare these 
redundancies in the sequel. The code over GF(2 bl ) is described 
since it is a bridging step to understand the one over GF(2). 

By Theorem [2] for the ^-polynomial 62(2;) — 1 + x + 
x 2 + ■ ■ ■ + x b2_1 , over GF(2 bl ), there exists an irreducible 
polynomial p 2 (x) of degree m 2 = r 2 — b 2 + 1 such that 
e 2 (x)p 2 (x) generates an optimum b 2 -burst-correcting code C* 



of length n 2 



_ (2*1)" 



2 —1 

- — and redundancy r2. Let 



C* = {f(x) e GF(2 bl )[x] : e 2 (x)p 2 (a:)|/(a:), deg/(z) < n 2 }, 



and let 

C 2 = {/(x)GGF(2 bl )N 



p 2 (x)\f(x), deg/(x) < n 2 }. 



The code C 2 is also of length n 2 and has m 2 redundancy 
symbols over GF(2 bl )- We will show now that C 2 can serve 
as the horizontal component code, i.e., given that the cluster 
occurred up to a cyclic permutation is e = (eo, ei, . . . , e(, 2 _i), 
where a G GF(2 bl ), for < i < 6 2 - 1, it will be possible to 
determine the first column i, < i < n 2 — 62, of the cluster 
with C 2 . 

Lemma 3: Let e = (eo, ei, . . . , e& a _i) be a given cluster, 
up to a cyclic permutation, which occurred in a transmitted 
codeword and found by the vertical component code. Then, the 
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(c) 

Fig. 2. Figure |2(a)| describes a cluster-error in a (2 X 3)-burst-correcting 
code. The cluster found by the vertical code is demonstrated in figure [2(b)] 
and figure [2(c)] shows the possible clusters considered by the horizontal code. 



horizontal component code C 2 can determine the first column 
of the given burst. 

Proof: We have to prove that if the burst e, or a cyclic 
shift of e, occurred in two different codewords fi(x), f 2 (x) 
then two different words will be generated. Since C 2 is a linear 
code it is sufficient to prove that there is no codeword which is 
equal to the difference of two clusters which are cyclic shifts 
of the cluster e. Assume the first column of the cluster e is i, 



i.e., the cluster is 



io i e io + l 1 



> e io+f>2-l 



where i a 



i(mod b 2 ) and indices are taken modulo b 2 . The polynomial 
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representing this cluster is 

hi(x) = x eib3 (e io x io + e i0+ ix i0+1 + ■ ■ ■ + e b2 ^ix b2 - 1 + 
e x b2 + e x x b2+1 + ■■■ + e^x^"' 1 ), 

where £j = L^J- We can write hi(x) as 

hi(x) = x Ub2 (e + e x x + --- + e h ,^ x x b2 ~ x + e (x b2 - 1)+ 
ei (x b2+1 - x) + ■ ■ ■ + e io _ 1 {x b * +i °- 1 - x 10 ' 1 )) = 
x l ' b2 (e + eix + ■ ■ ■ + e b2 - lX b2 - 1 + 
(x b2 - l)(eo + e x x + ■ ■ ■ + e^-ia:* - 1 )) = 
x lM (e + eix + • • ■ + e fc2 _nz; f,2 ~ 1 ) 
+x^ b2 (x b2 - 1) (e + e x x -\ h e^.ix 10 " 1 ) . 

Assume the contrary, that the difference between two clusters 
which are cyclic shifts of the cluster e is a codeword. Assume 
that these two clusters start at columns i and j. Hence, the 
polynomial hi(x) — hj(x) is the codeword 

hi(x) ~ hj{x) = x Ub2 (e + e lX -\ h eb^-ix^ 1 ) + 

x Ub2 {x b2 - 1) (e + e x x + • • • + e^-i^ " 1 ) - 
x^ b2 (e + e x x + • • • + e b2 - 1 x b2 - 1 ) - 
x^ b2 (x b2 - 1) (e + e x x + ■ ■ ■ + e^x^'- 1 ) . 

hi(x) — hj(x) can be written as 

hi(x)-hj(x) = {x Ub2 -x l i b2 ) 

(e + e x x + • • • + e b2 _ia; b2_:L ) + 
(x b2 - l)(x iib * (e + e x x + ■ ■ ■ + e^x 10 ' 1 ) 
-x 1 ^ 2 (e + e x x + ■ ■ ■ + e^-!^' " 1 ) ) = 
x^ix^- 1 ^ - 1) (e + e x x + ■ ■ ■ + e^x^- 1 ) + 
(x b2 - l)(x e > b2 (e + e x x + --- + e^x 10 - 1 ) - 
x^ b2 (e + e x x + ■ ■ ■ + e ]Q - X x J °- 1 )). 

This last presentation of hi{x) — hj(x) implies that it is 
dividable by x b2 — 1 and hence also by e2(x) = x ^Z x ■ Since 
hi(x) — hj(x) is a codeword in C2 it follows thatp2(a;)|^-i(a;) — 
hj(x). Since also e2(x)\hi(x) — hj(x), P2(x) is an irreducible 
polynomial, and its degree is greater than 62 it follows that 
e2(x)p2(x)\hi(x) — hj(x). Therefore, hi(x) — hj(x) is also a 
codeword of C*, a contradiction since C* can correct any burst 
of length 62- ■ 

A code C that can find the first column of a burst e = 
(eo, ei, . . . , eb_i) given up to a cyclic shift will be called a b- 
burst-locator code. Thus, by Lemma|3] C2 is a &2 -burst-locator 
code. Based on the constructions of [3], [8] and Lemma |3]we 
can construct an n x x ri2 two-dimensional (b x x ^-cluster- 
correcting code with small excess redundancy. 

Let C x be an optimum &i-burst-correcting code, over 

h (2 b 2\ r l— b l + 1 _ X 

GF(2 2 ), of length n x = - — 2 b 2-i redundancy r x . 

Let C2 be a 62 -burst-locator code, over GF(2 bl )> of length 
n 2 = 2 b i-i — anc ^ redundancy r7i 2 = r 2 — o 2 + 1. We can 
give a construction in which each codeword of size n x x 712 
has r x b 2 + rn^i + 1 = ri°g2( n i ri 2)l + b x b 2 + b x redundancy 
bits. The redundancy bits are partitioned into three subsets: 



• r x b 2 redundancy bits are located in the upper right corner 
of the array and are computed from the complete vertical 
codeword as done in Construction BBZS. 

> m2&i redundancy bits which are computed from the 
complete horizontal codeword as done in Construction 
BBZS. These redundancy bits are spread in the array in 
a way that they will fulfill the following requirement. 
If a redundancy bit is erroneous it will be possible to 
determine its row (note, that if the vertical code finds only 
one row where errors occurred, there are 61 different sets 
of 61 rows in which the burst occurred). Hence, in 2b x — 1 
consecutive rows there can be at most one redundancy bit. 
This requirement implies also that in any (61 x 6 2 )-cluster 
we have at most one redundancy bit. 

« One redundancy bit which is a parity of all redundancy 
bits of the second subset. Its role is to determine whether 
this bit or a redundancy bit from the second subset is 
erroneous. 



The encoding is done similarly to Construction BBZS 
with two exceptions. When we compute the elements of the 
horizontal component code the r x b2 redundancy bits of the first 
subset are not assumed to be zeroes as in Construction BBZS, 
but have the values which were computed by the previous 
steps of the encoding procedure. The second exception is the 
extra computation of the redundancy bit of the third subset, 
which is taken as an even parity bit of all the redundancy bits 
of the second subset. 

The decoding is done similarly to Construction BBZS with 
the following exceptions. In Construction BBZS, if redun- 
dancy bits are erroneous then they will be recovered by the cor- 
responding component code. The redundancy bits in the right 
upper corner are recovered by the vertical component code. 
They will be recovered by this code also in our construction. 
The redundancy bits in the left lower corner, of Construction 
BBZS, are recovered by the horizontal component code in 
Construction BBZS. Since, in our construction we don't use 
a burst-correcting code as an horizontal component code we 
cannot use this code to recover the related redundancy bits of 
the second subset. Each b x x 62 sub-codeword can contain at 
most one redundancy bit from the second or third subset. By 
summing all these redundancy bits we will know if one of 
them is erroneous. Also, these redundancy bits are spaced in 
a way that if we know in which rows errors occurred then we 
will know which one of these bits is in error. Once we will 
find this erroneous bit we will know the shape of the burst 
and the horizontal burst-locator code will find the column in 
which the cluster started. If only a redundancy bit from the 
second or third subset is in error then the vertical code will not 
find erroneous bits. In this case the sum of these bits is odd. 
If a bit from the second subset is in error then the horizontal 
burst-locator code will correct this error since this code is 
generated by a primitive polynomial and hence it can also 
correct a single error. Otherwise, the horizontal code will not 
find an error, which implies that the redundancy bit of the 
third subset is erroneous. 
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B. Binary burst-locator code 

The redundancy of the construction is improved if we use as 
the horizontal burst-locator code a binary (&i&2)-burst-locator 
code C of length 2 m — 1, where m = r — b\b 2 + 1, bib 2 

odd and e 2 (x) = 1 + x H H a; bl&2_1 . This is done simply 

by taking the 61 parity symbols which are computed for each 
column as 61 consecutive symbols in C 2 instead of an element 
in GF(2 bl )- 

Each codeword {cy } of size n\ x n 2 has r\b 2 + m + 1 
redundancy bits in the following positions: 

• ■ < i < n - l,n 2 - &2 < j < n 2 - 1}- 
These redundancy bits are computed from the vertical 
component code. 

. {(i(26i - 1) +j(2h - l)b x ,j) : < i < b x - 1, < 
j, (j + l)&i + i + 1 < m}. These redundancy bits are 
computed from the horizontal component code. 

• {(ni — 1, 112 — 1)}. This redundancy bit is an even parity 
bits for the redundancy bits of the second subset. 

Encoding: 

All the redundancy bits in a codeword are set initially to 
be zeroes. The vertical component code and the first set of 
redundancy bits are computed as in Construction BBZS, i.e., 
for each row i = r\, . . . , n\ — 1, 62 parity check bits are 
generated. pn, I = 0, 1, . . . , 62 — 1 is computed as 

PU = ^2 Cij. (2) 

j=e,e+b 2 ,e+2b 2 ,...,j<n2 

The parity bits pu, £ = 0, 1, . . . , 62 — 1 generate afterward 
a symbol p. — (p i0 jPii, ■ • ■ ,Pib 2 -i) from the extension field 
GF(2 &2 ). The symbols p., i = 7*1, . . . , n\ — 1 are considered 
as the information symbols of the code C\. By the encoding 
procedure for C\ we obtain r\ redundancy symbols p., i = 
0, ...,ri — 1, and the r\b 2 upper right corner redundancy 
bits of the array are computed in a way that (f2]i holds for 
i = Q, 1, . . . , n — 1 and £ = 0, 1, . . . , 62 — 1. We now turn 
to the encoding procedure of the horizontal component code. 
During this process the r\b 2 redundancy bits of the first subset 
will have the values which were just computed (as said before 
this is different from Construction BBZS, in which they were 
assumed to be zeroes). The second subset of redundancy bits 
spans over |"^] consecutive columns. In each column, with a 
possible exception of the last one, there are b\ redundancy bits 
in b\ positions (rows) which cover all the b\ distinct residues 
modulo b\. We compute rii — m information symbols of the 
horizontal component codeword as in Construction BBZS. 
The remaining m symbols are the redundancy symbols of the 
horizontal (foi&2)-burst-locator code and they are computed 
from the ri2 — m information symbols. The only redundancy 
bit of the third subset is the binary sum of the computed 
redundancy bits from the second subset. 
Decoding: 

The decoding is done similarly to the one of Construction 
BBZS. First, each row generates &2 parity bits by using (O 
(this includes also the first 7*1 rows, but the m + 1 redundancy 
bits of the second and the third subsets are assumed to be 
zeroes). These 62 bits are considered as a symbol in GF(2 &2 ) 
and hence a word of length m over GF(2 f ' 2 ) is generated. 



Now, we use the decoding procedure of the vertical 61-burst- 
correcting code to correct a burst of length b\. If such a 
burst occurred it "almost" determines the rows in which errors 
occurred and also the shape of the cluster up to horizontal 
cyclic permutation. We say "almost" since the vertical code 
does not find erroneous redundancy bits from the second and 
the third subsets. These bits are spaced in a way that at most 
one such bit is in error. We sum these m + 1 bits and if the 
result is not zero then one of these bits is erroneous. If this 
is the case then from rows of the cluster-error, discovered by 
the vertical code, we will know the exact row of this bit. If 
the vertical code didn't find any burst and a redundancy bit 
from either the second subset or the third subset is erroneous 
then we have exactly a single error in the array. Now, since the 
(&i&2)-burst-locator code is also a single-error-correcting code 
(binary Hamming code) this single error can be corrected. If 
there are more errors in the cluster then we continue by either 
correcting the redundancy bit of the last two subsets (and the 
corresponding bit in the horizontal codeword) or adding this 
erroneous redundancy bit to the shape of the burst. In either 
cases the horizontal (fri&2)-burst-locator code will discover 
the first column in which the cluster occurred, and hence the 
pattern discovered by the vertical component code enables us 
to correct the errors. 

Remark 1: The parity bits of the second set can be chosen 
in other ways as long as they form a set of redundancy 
symbols for the burst-locator code, e.g., they don't have to 
be in consecutive columns. Such choices can result in other 
array sizes. 

Remark 2: A natural question is to ask why not to use 
also a binary vertical component code? The answer is that 
we can. The main advantage will be that we will have more 
flexibility in the parameters of our two-dimensional array. The 
disadvantage is that the excess redundancy will be increased 
by [loga^l- 

The consequence of this construction is the following theo- 
rem (the computational part of the proof will be given in the 
next subsection). 

Theorem 3: The given construction produces a (bi x 62)- 

cluster-correcting code of size n\ x 712 = - 1 — 2t , 2 _ l x 

[- ^ — — with redundancy ri&2 + r — O1&2 + 2 < 

[log 2 (nin 2 )l +hb 2 + r io g 2 & il ■ 

The construction can be applied whenever &1&2 is odd integer, 
and there exists 

• An optimum b\ -burst-correcting code, over GF(2 b2 ), with 
redundancy ri and length i — 2 b 2-i ■ 

• A binary (61 62) -burst-locator code with redundancy r and 
length 2 r - blh2+1 - 1. 

C. Multidimensional arrays 

As said earlier, one of the advantages of our construction 
is that it can be generalized in relatively simple way to D 
dimensions, while the excess redundancy remains relatively 
small. Assume we want to construct a D-dimensional code of 
size ni x n 2 x • ■ • x rip which corrects a box-error of size 
bx x b 2 x • ■ ■ x bjj. Let B = \\f =l bi, where B is an odd 
integer. For the first dimension we use a component code of 
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length n\ over GF(2~) which corrects a burst error of size 61. 
In each of the other D — 1 dimensions we use a burst-locator 
code which locates the position of the burst and its cyclic 
permutation in the corresponding direction. In dimension i, 
2 < i < D, we use a bi -burst-locator code of length m over 

B_ 

GF(2 b i ). The code of the first dimension finds the position of 
the error in the first dimension and the shape of the error, with 
a possible cyclic shift in each of the other D — 1 dimensions. 
The code in dimension i, 2 < i < D, finds the location 
of the position where the cluster starts in dimension i. After 
each code discovers the position where the cluster starts in its 
dimension (note, that this can be done in parallel), we have 
the corresponding cyclic shift in each dimension of the box- 
error found by the first code. Hence, we can now form the 
actual burst which occurred and correct it. As before, we can 
use in dimension i, 2 < i < D, a binary B-burst-locator code. 
For the first dimension we can choose consecutive redundancy 
bits as the redundancy bits of the first subset. For each other 
dimension we will have to choose positions for the redundancy 
bits, which will fulfill the requirements for the redundancy bits 
of the second subset (only one redundancy bit will be needed 
as a parity bit for all the redundancy bits of this form for the 
burst-locator codes of all dimensions). 

Theorem 4: Assume b%, &2j • • • i &D are odd integers, B = 
YiiLi h an d the following codes exist: 

B 

« An optimum &i-burst-correcting code, over GF(2 b i ), with 

s 

redundancy n and length - . 

2 F T -1 

• For each i, 2 < i < D, a binary S-burst-locator code 
with redundancy r.- L and length 2 m ' — 1 = 2 Ti ~ B+1 — 1. 

Then, there exists a (bi x 6 2 x • • ■ x &d) -burst correcting code 



of size iii x Ti2 x 



I2 b i ) r i- b i+ 1 -l 

x n D = i '—g x 



and redundancy r\j- + ^f = 



1 < 



[log^nx • • • n D )] +B+ [log^^- 2 )] + 1. 

Proof: The existence of the code is implied by the 
proceeding description and we only have to compute its redun- 
dancy. For each j, 2 < j < D, rij = 2 ' ~ 1 

B 

(2^1" ) T1 _i> i+ 1 



> 



Therefore, mni ■ ■ ■ no > 



2"i 



n 



D 

j=2 



Now, 



■ (Ef=2 m 3 ■ 

log 2 (B D - 2 h) - 1. 



taking into account that for each j, 2 < j < D, rnj > B, and 
w.l.o.g. we can assume that for each i, 1 < i < D, bi > 1, 
we have that for each j, 2 < j < D, 2 m ibj > 2(D - l)B. It 
follows that log 2 (ni7i2 ' ' ' n D) > "x 1 — B+ (y^^ m,, — 1) 

Hence we have that the redundancy of the code is rij^ + 
Ef=2 m J + 1 < \\og 2 {n 1 ---n D )-\+B+\\og 2 {b 1 B D - 2 )\+\. 



When B is even we have to modify our method in order 
to obtain similar results. The modifications include binary 
component codes in all dimensions. Each one of the D — 1 
burst-locator codes locates the position of a cyclic burst of 
size B + 1. This modification is described in the next section. 



IV. Colorings for Error-Correction 

The constructions presented in Section [ill] are best applied 
when the volume of the box error is an odd integer. The reason 
is that Lemma Q] is true only when b is an odd integer. Hence, 
if the volume of the box-error is an even integer then the 
construction of Section|IlI]has to be used in a slightly different 
way. We have to apply the construction for correcting a box- 
error which has odd volume and contains the "real" box-error. 
The price will be an increase in the excess redundancy. In 
this section we offer a novel method which will be useful to 
correct a box-error with even volume and also for correcting 
other types of cluster-errors. The excess redundancy will be 
similar to the one of the constructions in Section [Til] 

A. The coloring method 

The constructions with binary component codes use D com- 
ponents codes from which the first one is a burst-correcting 
code and the other D — 1 codes are burst-locator codes. 
Position k in component code s is the binary summation 
of certain positions in the array, which were defined with 
correspondence to some modulo value related to s, k, and the 
D indices which define the position in the array. We generalize 
this idea to handle more complicated cluster-errors to a method 
which will be called the coloring method. A codeword is a D- 
dimensional array A (not necessarily a D-dimensional box). 
We want to correct any cluster-error with a given shape whose 
volume is B. 

Again, we use D binary component codes to correct the 
cluster-error. The first code is a (B+Si) -burst-correcting code, 
5i > 0. The s-th component code, 2 < s < D, is a (B + 5 S )- 
burst-locator code, S s > S\. We further use D different 
colorings of the D-dimensional array. To each position of 
A we assign a color for each one of the D colorings. Each 
coloring will be associated with a different binary component 
code. For a given coloring and the corresponding component 
code, position k in the component code is the binary sum of 
all bits which are colored with color k. As we want to correct 
a cluster-error of a certain shape in the array we want that the 
colorings will satisfy a few properties: 

• (p.l) For the s-th coloring, for each s, 1 < s < D, 
the colors inside a burst of the given shape are distinct 
integers and the difference between the largest integer 
and the smallest one is at most B + 5 S — 1. 

• (p.2) Given D colorings, a color v s for the s-th coloring, 

1 < s < D, There is at most one position in the array 
which is colored with the colors {v\, ^2, ■ ■ ■ , vd). 

• (p.3) Any two positions which are colored with the same 
color by the first coloring, have colors which differ by 
a multiple of B + 5 S by the s-th coloring, for each s, 

2 < s < D. 

Finally, we have to choose redundancy bits in the array, in 
a similar way to the method used in Section [ill] 

Theorem 5: Assume that there exists a (B + <5i ^burst- 
correcting code of length n\ and for each s, 2 < s < D, there 
exists a (B + S s ) -burst-locator code of length n s . Assume 
further that there exist D colorings which satisfy properties 
(p.l), (p.2), and (p.3), such that the s-th coloring assigns colors 
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between 1 to n s to the D-dimensional array. Then the coloring 
method implies the existence of a £>-cluster-correcting code 
for a D-dimensional array A and a cluster with a given shape 
and volume B. 

Proof: The proof is straightforward from the description. 
We just note, that by property (p. 2) there is no ambiguity in 
the erroneous positions. By property (p.l), for each s, 1 < s < 
D, the erroneous positions affect at most B + S s consecutive 
positions in the codeword of the s-th component code. Finally, 
since two positions in the array which are assigned the same 
color by the first coloring, have been assigned by the s-th 
coloring colors which differ by a multiple of B + S s (see 
property (p. 3)), it follows that the possible bursts in the s-th 
(B + <5 s )-burst-locator code are cyclic shifts of a burst with 
length B + 5 S . ■ 
It will be more convenient if each coloring is a linear 
function of the coordinate indices, i.e., given a position 
(£i,i2, • • • its color for the s-th coloring will be defined 
by 

D 

E aS i l J 

where a? is a constant integer which depends on the coloring 
s and the shape of the £>-dimensional cluster. Such a coloring 
will be called a linear coloring. With a linear coloring we 
associate a coloring matrix Arj, where (An) s ,j — a j- It 
is easy to verify that property (p.2), is fulfilled for a linear 
coloring if and only if the coloring matrix is an invertible 
matrix. 

One can verify that the coloring method is a generalization 
of the method described in Section [HI] To observe this we 
should assign a color to position (£i, i%, . . . by the s-th 
coloring, 1 < s < D, in a slightly different way than the 
assignment in the next subsection. If position k in component 
code s is the binary sum of a certain set S of positions in the 
array, then all positions of S are colored with color k by the 
s-th coloring. We leave the exact definitions of the colorings 
as an exercise for the reader. 



Each coloring corresponds to one component code. Codeword 
of component code s, C s , 1 < s < D is defined according 
to coloring s. Position k in the codeword is the sum modulo 
2 of the values in positions colored with color k by the s-th 
coloring. We will prove that these D codes satisfy properties 
(p.l) through (p.3). 

Lemma 4: If the D-dimensional code has a box error of 
size bi x 6 2 x • • • x brj then each one of the D component 
codewords has a burst whose length is at most B. 

Proof: A cluster occurred in the array is located inside a 
multidimensional box of the form {(£* + i\, £| +i2, ■ ■ ■ , ijj + 
id) '■ < ij < bj — 1, 1 < j < D}, for a fixed position 
(£*, »2, . . . The smallest color, lc s , of an erroneous 

position in C s , 1 < s < D is 



s-l 



1) 



D R 

E 



3=s 



which is the color of position + b% — 1, . . . , i* s _\ + 6 s _i — 
l,i*,...,i D ) and the largest color, hc s , of an erroneous 
position in C s , 1 < s < D is 

^=E-^- 1 |^£*+E^+^- 1 )' 

, J-'s — 1 J-Js— 1 

which is the color of position (£*,..., £*_ x , £* + b s — 
1 J, + & d — !)■ Now, we compute the difference hc s — ic s ■ 

h Cs ~ ic. - E B J-i4^^ 1} + E f^fe - !) 

-. £>s — 1 -°s — 1 



s-1 £> 

E>-i& - 1)) + -— (E^-ife - 1) 



Bp 

Bj, 



-(B s -!-B ) 



B. 



s-l 



B D - S s _i 



B. Multidimensional box-errors 

To demonstrate how the coloring method works we will first 
show how it is used to correct multidimensional box-errors, 
where the volume of the box error is an even integer. 

Assume we want to construct an rii x ra 2 x ■ • ■ x nn 
D-dimensional x t 2 x • ■ ■ x &£>)-cluster-correcting code, 
where B — YliLi h is an even integer. We will use D binary 
component codes. One component code will be able to correct 
a burst of length B and D — 1 component codes will be able 
to locate the position of a burst, whose length is B + 1, given 
by a cyclic shift. Let 



= B D -^ + J^-l = B-l 

Therefore, the length of a burst in each component code is at 
most B. ■ 
Lemma 5: For each one of the D colorings, the B colors 
in each D-dimensional box of size b\ x 62 X • ■ ■ x bu in the 
array are all distinct. 

Proof: Assume the contrary, that there exist two different 
positions . . . , io), (ti, • ■ • , to) located inside a box of size 
b± x 62 x ■ ■ • x bo in the array whose s-th color is identical. 
Therefore, by definition, 



B 4 



bjBj-i, where 1 < j < D, and Bo = 1 



For each entry («x, £2, • • • > £u) in the array we assign D 
colors. The s-th color, 1 < s < D is defined by 

s-l „ D 



a U 2 -^D - E ^J-: 1 



B s _ 



B, 



B s -i 



3=1 



B 



B7-1 



E /; .- 1 77-^7 '• + E 77 



1 r *i + 2^ r *J 

£>s-i £>s-i 

3=1 J = s 



(3) 
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which implies 



E / ''' 

3=1 



J=«+l 



-Bs-i 



5/ 



6 S [ij-tj) 



Bj-i 



(4) 



Since these two positions are located inside a box of size 
b\ x • • • x bp it follows that < |f s — i a | < b s — 1 and hence 
f s — i s = 0. We continue with and by induction we prove 
similarly that tk~ik — for each k, s+1 < k < D. Therefore 
by we have 



5 



which implies 



2- J B i _ 1 (i J -t i )=0 

3=1 



Now, we will show by induction that for each j, 1 < j < s—1, 
tj — ij = 0. For j = 1, we have 



S— 1 S — 1 

tx-h =J2Bj-i(ij ~tj) = 6i(E -^fe - 



i=2 



J=2 



and since < |ii — ii| < b\ — 1 we have tx — ix = 0. We 
continue in similar way and obtain for each j, 1 < j < s — 1, 



0. Therefore, for each 1 < s < D, we have t s 



Thus, for each one of the D colorings, the B colors in each 
D-dimensional box of size ^ x t 2 x • ■ ■ x ^ in the array are 
all distinct. ■ 

Lemma 6: Any two positions which are colored with the 
same color by the first coloring, have colors which differ by 
a multiple of B + 1 by the the s-th coloring, for each s, 2 < 
s < D. 

Proof: Assuming the fc-th position in the codeword of C\ 
is erroneous. This error results from an array error in position 
(ii,«2, ■■■,*£>) such that a\ ii2 ^ 



i D = Ej=i Bj-xij = k, and 
hence ix = k — $^.-_ 2 Bj-xij- The possible error locations in 



C s , 2 < s < D are of the form 

s-X 

a iii 2 ---iD ~ 

i=i 

B D 



E IS > 7^' 



D R 



B„ 



B s -x 

3 = s 



D 



(k-j2 B j-iij)+J2- B ^ 



B D 



D V" / i'.;/ i / . — j — i R l 3 

E>s — X £>s — X 

J= 2 j=2 

D R 

^ B 



3=s 



s-X 



Br 



B s -x 
D R 

E^ 



3=s 



B s . 

Bp 
B s _x 

Bp 



J=2 ^ s - 1 



3=s 



s-l g 

E^-i^:^ 

3=2 ° S - X 



B 



i-i> 



#s-l 



s-l 



^ p k is a constant and therefore, two positions which have 
the same color by the first coloring, have colors which differ 
in a multiple of B + 1 by the s-th coloring, 2 < s < D. ■ 
Lemma 7: Given a position c s in the s-th component code, 
1 < s < D, the set of equations 



al li2 ... iD ,l<s<D 



(5) 



has exactly one solution for (ix, ■ ■ ■ , ip). 

Proof: We will prove that the coloring matrix is invertible 
by proving that its determinant is not equal to zero. The (s, j) 
entry of the coloring matrix Ap is given by 




for j < s 
for j > s 



For the basis of the induction, A2 



and hence 



We will prove by induction on D, D > 2, that \Ap\ — (1 
S) D_1 , where B is the volume of the box-error. 

1 bx 

-b 2 1 

IA2I = I + 6162 = I + 6162, where the two-dimensional cluster 
has size 61 x 62. 

For the induction hypothesis we assume that the determinant 
of every coloring matrix of size (D — 1) X (D — 1) is given 
by |Ad_i| = (1 + B') D ~ 2 , where B' is the volume of the 
corresponding (D — l)-dimensional box-error. 
For the induction step let Ap be a coloring matrix of size 
D x D. The determinant of Ap is given by 

p 

\A D \ = Y J {-W +1 Bj-i\A D [l,j]l 
3=1 

where Ad[1,j] is the matrix obtained from Ap by deleting 
row 1 and column j. For 2 < s < D — 1 the s-th row of 

j4.£)[1, s] is given by: 

-(—BpB , —BpBx, ■ ■ ■ , —BpB s -2-,B s . . . . , Bp-x), 



B 



s-l 
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and the (s + l)-th row is given by: 
1 



B 



■(—BdB , —BjjBi, . . . , —Be>B s -2, B s , . . . , Bd-\). 



These two rows are linearly dependent and therefore, 
|(4d)i s | = 0for2<s<D-l and we have 

\A D \ = B \(A D )u\ + (-l) D+1 B D ^\(A D ) 1D \. (6) 

The matrix (Ad)ii is also a coloring matrix with respect 
to coloring related to the box-error x x • • • x bo-\ X 
(bob\) and according to the induction assumption its deter- 
minant is given by (Ad)h = (1 + B) D ~ 2 . Let (A' D )i D 
be the matrix constructed from (Ad)id by dividing each 
element of the last row of {Ad)xd by — bp and shifting 
this row cyclically to be the first row. Clearly, |(Ad)i.d| = 
— bo{— l) D ~ 2 \(A' D )iE,\. The matrix (A' d )id is also acoloring 
matrix with respect to coloring related to the error box 
bi x b 2 x ••• x b D - 2 x (b D -ib D ). Hence, \(A d )id\ = 
-b D (-l) D - 2 (l + B) D - 2 and from © we have \A D \ = 
(1 + B) D - 2 + (-l) D+1 B D ^(-b D )(~l) D - 2 (l + B) D - 2 = 
(1 + B) ^ 1 . Thus, Ad is invertible and (0 has a unique 
solution. ■ 
The encoding procedure is quite straight forward. First we 
have to choose three sets [or D + 1 sets] of positions for the 
redundancy bits as in the previous constructions. Position k in 
the codeword of the s-th component code is the binary sum of 
all positions in the array colored with k by the s-th coloring. 
In the decoding procedure, the first component code provides 
a list of erroneous positions. If k is the erroneous location of 
the first code, the difference between equivalent errors in C s , 
2 < s < D is a multiple of B + 1 and by Lemma [6] the value 
— ® D k is the residue modulo B + 1 of the error location 

i>s-l 

in the codeword of C s . Therefore, for the codeword of C s the 
burst is known up to a cyclic permutation of length B + 1. 
C s is a burst-locator code which can locate where such burst 
started. Hence, each burst-locator code locates the position in 
which the burst has started with respect to the corresponding 
dimension. Therefore, the location of the erroneous position k 
in the first code is known in each of the D — 1 burst-locator 
codes. Hence, we can partition the positions of the errors in 
the D component codewords into D-tuples. Each such tuple 
is of the form (c\, c 2 , . . . , cd), where c s is the error location 
in the codeword of C s , 1 < s < D. By Lemma [7] such a D- 
tuple corresponds to exactly one position in the _D-dimensional 
array, which is an erroneous position. 

To summarize, by Lemmas |4] and [5] the D colorings satisfy 
(p.l), by Lemma [6] they satisfy (p. 2), and by Lemma [7] they 
satisfy (p. 3). Hence, by Theorem [5] the code is a (&i x b 2 x 
• • x bo) -cluster-correcting code. The redundancy of the code 
is slightly larger than the one for odd B (see Theorem H}. We 
omit the tedious proof. 

V. Lee Spheres Cluster Errors 

An error event at a position (ii, i 2 , ■ ■ ■ , id) can affect other 
positions around it. If we assume that the error event can 
be spread up to radius R, then any position (t\,t 2 , ■ ■ ■ , ir>) 
such that \tg — ie\ < R might be erroneous. The set of 

positions {(ti, t 2 , ■ ■ ■ , i_o) : Y^i=i \U ~ — R} forms a 



D-dimensional Lee sphere with radius R [6], [10]. Another 
important observation is that any arbitrary cluster-error of 
size b is located inside a Lee sphere with radius [|J- We 
wish to supply a method to correct all errors which might 
occur in such a given Lee sphere. The size of such sphere is 



E 



mm{D,R} r,j (D 



2 J [10] and each position can be erroneous. 



We can think about three different methods to perform the task. 

• We can use any method which is able to correct all errors 

occurred in a (2i? + 1) x ■ • • x (2R + 1) Z?-dimensional 
< ' 

D times 

cube since a Lee sphere with radius R is located inside 
such a cube. The main disadvantage of this method is 
that the code can correct much more errors than the 
ones needed in its goal. For example, the size of the 
Z?-dimensional cube is (2R + 1) D , while the size of 

the D-dimensional Lee sphere with radius R is ■ L -pj h 

OiR ' 1 ) when D is fixed and R — > oo. 
> We can transform the space into another space in a way 
that a Lee sphere will be transformed to a shape located 
inside another error-shape with a small volume which we 
know how to correct efficiently. 

• We can find a coloring and perform correction of errors 
as explained in Section [TV] 

In this section we will explain how we apply the last two 
methods. 

A. Space transformation 

The idea we are going to use is to transform the space into 
another space in such a way that a Lee sphere in one space 
will be transformed into a shape located inside a box in the 
second space. After that we will be able to use the encoding 
and decoding introduced in Section [Til] We start with the two- 
dimensional transformation. 

Lemma 8: Let M, M* be infinite two-dimensional arrays 
and let T be the transformation from M into M * defined by 
T (ii,i 2 ) = ( P 1 ^"' 2 ] , i 2 — ii). Then a Lee sphere with radius 
R in the array M is located after the transformation T inside 
a rectangle of size (R + 1) x (2R + 1) in M*. 

Proof: A Lee sphere of radius R with center at 
(ii,i 2 ) in the array M includes the set of positions 
B R (ii,i 2 ) = {(h,t 2 ) : |*i - ii| + \t 2 - *a| < R} = 



{(ii+R e ,i 2 + Rj) 



Rt. 



\Rj\ < R}. B R (i u i 2 ) is 



transformed by T into the set of positions B R (ii,i 2 ) 



{( 



h+Rt+12+Rj 
2 



.12 



of M*. Denote, i\ 
we have that 

{(i*x + t u q + t 2 



B*r' 



l-i? 



j ^i Re 



«2 - n 



\Rj\ < R\ 

R, and 



2 I ' "2 

ii,i 2 ) is located inside the rectangle 
: < h < R,0 < t 2 < 2R}. M 
The transformation T transforms a parallelogram into a 
rectangle (see Fig. [3] and @]l and hence we will need some 
adjustment in our encoding and decoding procedures if we 
want to correct Lee sphere clusters in a rectangular array rather 
than a parallelogram. 

The construction is generalized to D dimensions. 
The transformation T will work between two D- 
dimensional arrays. For each entry (ii,i 2 , . . . ,Id) in the D- 
dimensional array M, T (ii,i 2 , . . . ,id) = (ij:i 2 , • ■ ■ i *£>)> 
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Fig. 3. Positions in the parallelogram 
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Fig. 4. Transformation to rectangle 



where ? 

(-i) 



= m 



= r- 



n+ - '2+»3 

2 



for 1 < j < D - 1 and 



1 — I ~ 2~ I ' '2 

J + 1 n + (-l) 3 ^2-l + 1 

2 

Jjj = Y>® =l (—l) D ~Hj. We invoke first the two-dimensional 
transformation T on the first two coordinates, then on the 
second and the third coordinates and so on. Similarly to the 
proof of Lemma [8] we can prove the following lemma. 

Lemma 9: Let M be an infinite D-dimensional array. Then 
a Lee sphere with radius R in the array M is located after 
the transformation T inside a D-dimensional box of size 
(R + 1) x (R + 1) x • • • x (R + 1) x(2R + 1) in M*. 

s / 

73—1 times 

As a consequence of this transformation we can use the 
constructions of Sections [Til] and [IV] for correction of D- 
dimensional box-error. If we assume that our codewords are 
£>-dimensional arrays rather than D-dimensional parallelo- 
grams we can use an array located inside the parallelogram. 
The redundancy in this case will be slightly increased. 

B. Tiling and coloring 

When the error shape is a two-dimensional Lee sphere with 
radius R we can provide a code with a better redundancy 
than the code constructed by using the two-dimensional space 
transformation. We will use the coloring method of SectionHVl 
Indeed, this construction is a good example for efficient uses 



of the coloring method. We choose two appropriate colorings 
and , J , 2 . For a given position (ii,i 2 ) m the n\ x n 2 
two-dimensional array let ^x{ii,i 2 ) = {R + l)*i + -R*2 and 
*2(*i, 12) = -Rh + {R+ l)i2 (see Fig. B- 

Lemma 10: If \&i(ii,«2) = ^i(h,t 2 ) then ^2{h,h) - 
^2(h,t 2 ) is a multiple of 2R 2 + 2R + 1. 

Proof: *i(ii,i 2 ) = ^i(ti,t 2 ) implies that (i? + l)ii + 
Ri 2 = (R + l)t 1 + Rt 2 . Hence, R(i 2 -t 2 ) = {R+ l)(ti -ii), 
i.e., (R+l)(i 2 -t 2 ) = £±p^(t 1 -i 1 ). 

Now, * 2 (ji,i 2 )-*2(ii,t 2 ) = -Rii 
(R + l)t 2 ) = R(h - h) + (R+ l)(i 2 



-(R+l)i 2 -(- 
- 1 2 ) = R(h - 



Rh + 
h) + 



R 

Since R and 2R 2 + 2R + 1 are relatively primes it follows 
that R divides t\ — i\ and hence vE , 2(ii,* 2 ) — ^2(^1,^2) is a 
multiple of 2R 2 + 2R+1. ■ 

The effect of these colorings of the two-dimensional arrays 
is best seen if we consider a tiling of the two-dimensional 
space with Lee spheres with radius R. Such tiling is well 
known and given in [6], [10]. In this tiling each Lee sphere 
belongs to two diagonal strips. For the first coloring all relative 
positions of the Lee spheres in the same diagonal strip have 
the same number; and in the other direction they are congruent 
modulo b* = 2R 2 + 2R + 1 which is the size of a sphere (see 
Fig.0. 

Lemma 11: For any given Lee sphere with radius R in a 
two-dimensional array, the color numbers assigned by ^>i, i = 
1,2, to the positions of the Lee sphere are 2R 2 + 2R + 1 
consecutive integers. 

Proof: Given a Lee sphere centered at the point (ii,i 2 ) it 
is readily verified that smallest color number assigned by ^1 
to a point in the Lee sphere is — R, i 2 ) and the largest 

color number is + i 2 ). Now, ^(ii + R, i 2 ) — 4 r i(*i — 
Ri 2 ) = (R + l)(i 1 +R) + Ri 2 -((R+l)(i 1 -R) + Ri 2 ) = 
2R 2 + 2R. Hence, to complete the proof we have to show that 
all the color numbers assigned by ^1 to the positions inside 
a Lee sphere are distinct. 

Let (ii,i 2 ) and (tx,t 2 ) be a pair of points for which 
^liii)^) — ^1 (ii , ^2 )■ By the proof of Lemma [TOl we have 
that R(i 2 - t 2 ) = (R + - Hence, i 2 - t 2 = 

-^(ii - ii) and t 1 -i 1 = j^ifa - h)- Since R and R+ 1 
are relatively primes it follows that R divides t\ — i\ and 
R + 1 divides « 2 — < 2 . This implies that R < \ti — ii\ and 
R+ 1 < \t2 - t 2 \. Therefore, 2R + 1 < \t% - ii\ + \i 2 - te\ 
and the two points (ii,i 2 ) and (ti,t 2 ) cannot be contained in 
the same Lee sphere. 

Thus, all color numbers assigned by to the positions 
contained in a Lee sphere are 2R 2 + 2R + 1 consecutive 
integers. The same proof holds for 4' 2 . ■ 

Lemma 12: The coloring matrix defined by "fi and 4' 2 is 
an invertible matrix. 

Proof: The coloring matrix defined by ^1 and ^ 2 is 

R+l R 
-R R+l 

which is clearly an invertible matrix. ■ 
By Lemmas 1101 [TT1 .and [T2l we have that (p.l), (p. 2), and 
(p. 3), are satisfied respectively. Hence, by Theorem[5]the code 
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constructed is capable to correct a Lee-sphere error with radius 
R and size b* . 

We choose r+m+l appropriate redundancy bits similarly to 
the constructions in Section Hill We use two component codes. 
One code is a b* -burst-correcting code of length 2 r ~ b +1 — 1, 
where 2 r ~ b +1 is the least power of 2 greater than (R+l)(rix— 
1) + R(n 2 — 1) + 1 which is the number of colors needed to 
color the array with the first coloring. The second component 
code is a b* -burst-locator code of length 2 m — 1 (such a code 
exists since b* = 2R 2 + 2R + 1 is an odd integer). We take 
2 m to be the least power of 2 greater than (R + l)(ri2 — 
1) + R(ni — 1) + 1 which is the number of colors needed 
to color the array with the second coloring (for the simplicity 
of the computations we will take 2 r ~ b +1 > (R + l)ni + 
Rn 2 > 2 r - b ' and 2 rn > (R + l)n 2 + Rni > 2™" 1 ). The 
redundancy of the n\ x n 2 two-dimensional code is at most 
[log 2 (raira 2 )l + b* + r21og 2 (2i? + 1)] + 2. If n x = n 2 = n 
then the redundancy of the two-dimensional code is at most 
|~log 2 n 2 ] + b* + |~21og 2 (2i? + 1)] and if each codeword is a 
rhombus of size n 2 then the redundancy is at most |~log 2 n 2 ] + 

&*+riog 2 &*i. 

C. Multidimensional codes 

A similar idea can be used for L>-dimensional code of size 
ni x ri2 x • • • x ud correcting Lee sphere error with radius 
one. We use D different colorings of the array. For a position 
(i x , i 2 , . . . , ip), < ie < n — 1, 1 < £ < D, the s-th 
coloring, 1 < s < D, assigns the color ^f =1 £ ■ i s +t-i, where 
indices are residues modulo D between 1 and D (note that 
each color in each coloring forms a perfect code [10] when 
we consider the coloring in lP and not just in the array). 
Again, for each coloring, all the color numbers located in a 
Lee sphere of the Z?-dimensional array are 2D + 1 consecutive 
integers. Bit k of the s-th component code is the binary sum 
of all the bits colored with the integer k, by the s-th coloring, 
in the Z?-dimensional codeword. It is easy to prove that the 
coloring matrix is invertible. But, property (p. 3) does not 
hold in all dimensions. Therefore, if property (p. 3) does not 
hold for the s-th coloring we need to take a (2D + 1) -burst- 
correcting code instead of (2D+ 1) -burst-locator code for the 
s-th dimension. The consequence will be a code with larger 
redundancy. Clearly, there are many appropriate colorings for 
each dimension. Hence, we can try to replace each coloring 
by one for which property (p. 3) holds. Of course we have 
to make sure that the coloring matrix will be invertible. The 
excess redundancy in this case is quadratic in D, compared to 
exponential in D if we use the transformation T and the code 
which corrects D-dimensional box-error. 

VI. Bursts with Limited Weight 

We now turn for a new kind of errors which are important in 
our final goal of correcting an arbitrary cluster-error. Assume 
we need to correct errors in a cluster of size b, where the 
number of erroneous positions is at most t. We wish to find 
one-dimensional and multidimensional codes which correct 
such errors. These codes have an obvious application as we 
can expect that in area that suffers from an event which caused 



errors, not all positions were affected. Another important 
observation is that any arbitrary cluster-error of size b is 
located inside a Lee sphere with radius [|J- Hence, this 
cluster-error can be corrected if we can correct a Lee sphere 
error with radius |_§J , where the number of erroneous positions 
is at most b. 

A. One-dimensional codes 

Throughout this subsection all codes are binary and & is an 
odd integer. 

Let Ci be a i-error-correcting code of length b, and redun- 
dancy r. Let Hi be its parity-check matrix of size r x b, 

Hi = [hi, h\, . . . , h\_i\. 

Let C 2 be a 6-burst-locator code, of length n = 2 m — 1, and 
redundancy m. Its parity-check matrix, H 2 , is of size m x 
(2 m - 1), 

H 2 = [h , hi, . . . , h 2 n ^i\. 

Based on these two matrices, we construct a new parity- 
check matrix H, 

h{ ■ 

hi ■ 

The j-th column, < j < n 



H = 



hi 

>4 



hi 



h) 

(n 



-r)(mod b) 



"n-1 

1, of the matrix H will be 



defined as the concatenation of the j(mod b)-th column of the 
matrix Hi and the j-th column of the matrix H 2 , 

Lemma 13: H is a parity-check matrix for a code C of 
length 2 m — 1, correcting every burst of length b with at most 
t erroneous positions. 

Proof: It is sufficient to prove that there are no two 
bursts of length b and weight at most t that have the same 
syndrome. Assume the contrary, that there are two words, of 
weight at most t, with bursts of length b, y 1 = (j/q, . . . , 
and y 2 — (j/q, . . . , y 2 _i), which have the same syndrome, i.e. 
H(y 1 ) T = H{y 2 ) T . This implies that 



n-l 



n-1 



(7) 



i=0 
n-l 



i=0 
n-l 



E»fa a = E*fo 9 > (8) 

i=Q i=Q 

where [j]e is the unique residue of j(mod £) between and 
£ — 1. Let %l, be the first nonzero bit in y 1 , y 2 , respectively. 
Therefore, (0 can be written as 

il+b-l il+b-1 
i—iQ i— Iq 

If we denote £ Y = b ,£ 2 = b then the last equation 
can be written in the following way: 

[jjjb-l 6-1 

yl'+b+M + yp+ihl = 

i=0 

yi 2 +b+i h l+ y 2 i 2 +* h i- 

i=0 <=[ig]t 



b-1 
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(a) (b) 
Fig. 5. (a) The coloring ^i, and (b) the coloring ^2. with R = 2. The left upper corner coordinate is (—6, 9). 



This equation implies that for the words of length 6 

zl = (vh+b* ■ ■ • ■ ■ ■ >yh+b-i) « 

2(2 2 2 2 \ 

we have Hi(z 1 ) T = Hi(z 2 ) T . The weight of these words is 
at most t and since Hi is a parity-check matrix of a t-error- 
correcting code we have that z 1 — z 2 . Therefore, the two 
different words y 1 and y 2 contain the same burst of length 6 up 
to a cyclic permutation. By Lemma[3]we have that H2(y 1 ) T ^ 
H2(y 2 ) T , contradicting ®. 

Thus, H is a parity-check matrix for a code of length 2 m — 1, 
correcting every burst of length & with at most t erroneous 
positions. ■ 
Remark: It is important to note that H is acting similarly to 
the constructions of previous sections. The first part which is a 
concatenation of several copies of Hi is a parity-check matrix 
of a 6-burst-correcting code with at most t erroneous positions. 
It finds the burst pattern up to a cyclic shift. The second part 
of TC, H2, is a parity-check matrix of a 6-burst-locator which 
can find the location of a burst of length b given up to a cyclic 
shift. Lemma Qj] gives a formal proof for these facts in terms 
of the generated syndromes. 

To summarize the parameters of the construction we need 
the parameters of t-error-correcting codes of length b. We 
can use BCH codes [13], [15] for this purpose. If I is the 
least integer such that b < 2 e — 1 then there exists a i-error- 
correcting BCH code of length 2 l — 1 and redundancy at most 
t£. By shortening we can obtain a t-error-correcting code of 
length b and redundancy t\log 2 b~\. Now, we can summarize 
our construction in this section. 

Theorem 6: The code C has length n = 2 m — 1 and it 
corrects every burst of length b with at most t erroneous 
positions. If b is an odd integer then the redundancy of the 



code is m + t\log 2 b~\ < |~log 2 n] + i[log 2 6], and if b is an odd 
integer then the redundancy of the code is m+t [log 2 (b+ 1 )] < 
[log 2 nl+iflog 2 (&+l)l. 

B. Multidimensional codes 

Now, we want to design a two-dimensional code of size ni x 
ri2 capable of correcting a (61 x 62) -cluster with weight at most 
t, where 6162 is an odd integer. We use a construction similar 
to the one used in previous sections. In this construction we 
have to use two binary component codes. The vertical one is 
a (61 62) -burst-correcting code, of length 72462 in which the 
weight of the burst is at most t. Such a code was constructed 
in the previous subsection. The horizontal component code 
is a (6i&2)-t>urst-locator code of length n^bi. They are used 
in the same manner as they are used in previous sections to 
correct a burst of size bi x 62. Since the vertical code can find 
a burst only if the weight of it is at most t, it follows that the 
two-dimensional code can handle bursts of size bi x 6 2 only 
if their weight is at most t. If 6162 is an odd integer than by 
Theorem [6] we can take a vertical code of length 71162, 2 m > 
72462 > 2 m_1 , with at most m+t [log 2 (6162)] redundancy bits. 
The horizontal code has length 2 r - blb2+1 - 1, 2 r - blb2+1 > 
71261 > 2 r ~ blb ' 2 , with r — 6162 + 1 redundancy bits. Therefore 
we have. 

Theorem 7: The redundancy of the m x code C which 
is capable to correct a (61 x 62)-cluster, 6162 an odd integer, 
with weight less or equal t, is at most [log 2 (77,1712)] + it + 
l)riog 2 (6 1 6 2 )l +3. 

The generalization for multidimensional codes is straight 
forward. If the size of the cluster is an even integer then we 
will use an appropriate coloring. 

VII. Correction of Arbitrary Bursts 

Finally, we want to design a code which corrects an arbitrary 
£>-dimensional cluster-error of size 6. If 6 is odd then the 
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cluster is located inside a Lee sphere with radius If b is 
even then either we consider it as a cluster-error of size 6+1 
or slightly modify the constructions for a small improvement 
on the efficiency. Modification based on coloring can be also 
obtained if 6 = 3 and D > 3. 

We now show how the codes of Section [VI] help to correct 
an arbitrary cluster-error of size b. We start again with two- 
dimensional codes. A cluster-error of size b is located inside a 
6x6 square. Therefore, we can use a code of size n 1 xn 2 which 
corrects a (6 x 6)-cluster with weight 6. The generalization for 
D-dimensional code is straight forward. 

An improvement in the excess redundancy is obtained if we 
consider a code which corrects a smaller shape, with limited 
weight, in which the 6-cluster is located. For simplicity we 
will consider only the case where 6 = 2R + 1 and the small 
shape is a Lee sphere with radius R. We will use the colorings 
and ^2 given in Section [VI and two component codes, the 
first one is a 6* -burst-correcting code, 6* = 2R 2 + 2R + 1, in 
which the weight of the burst is at most 6, and the second one 
is 6* -burst-locator code. Now, we apply the coloring method 
of Section [TyJ The redundancy computations are similar to 
the ones in subsection IV-BI and in Theorem [7] and we have the 
following result. 

Theorem 8: the redundancy of an n\ x n 2 code capable to 
correct a cluster of size 6 is at most |~log 2 (ni?i 2 )] + (6 + 
l)r21og 2 61+3. 

Generalization for £>-dimensional code is done by using 
the transformation T of Section [Vj Hence, there is some 
loose of efficiency, but the performance is still better than the 
performance of a code which corrects a D-dimensional cluster 
error whose shape is a Z?-dimensional box-error with limited 
weight 6. The redundancy computation is similar to the ones 
in previous sections. 

The next question of interest is a lower bound on the excess 
redundancy of a code which corrects an arbitrary cluster of 
size 6. A lower bound on the excess redundancy is log 2 iV£>(6), 
where Nrj (6) is the number of distinct patterns considered as 
D-dimensional clusters of size 6. Finding bounds on Nrj(b) 
is an interesting geometrical combinatorial problems of itself. 
A related question is to find number of distinct clusters with 
size 6, with no "holes", and exactly 6 erroneous positions. This 
problem is the same as finding the number of 6-polyominos. 
For D — 2 the known lower bound on their number is 
3.981037 b [5] and the known upper bound is 4.649551 fa [12]. 
Therefore, we have 

Theorem 9: The excess redundancy of a two-dimensional 
code, which is capable to correct an arbitrary cluster of size 
6, is at least 6 • log 2 3.981037 . 

Theorem[9]is a small improvement of the trivial lower bound 
(which is 6) on the excess redundancy. But, the gap between 
the orders of the lower bound 0(6) and the upper bound 0(6 • 
log 2 6) is still large. 

VIII. Representation with Parity-Check Matrices 

All codes which were discussed in the previous sections 
are using auxiliary linear codes for the computation of the 
redundancy bits in the codewords and to conduct the proper 



decoding. It is not difficult to see that all the multidimensional 
codes are linear, by noting that the bit by bit addition of two 
codewords is another codeword. In this section we will explain 
how to present similar codes with parity-check matrices. This 
will be done by considering all the component codes as binary 
codes. 

Remark: When the component code is a linear burst- 
correcting code of length n over GF(2 b ) we can consider it as 
a binary code of size n x 6, as we actually use it. We note 
that the bit by bit addition of two codewords of size n x 6 is 
also a codeword and hence the code is a binary linear code. 

In these new codes, which will be constructed, we won't 
need the redundancy bit of the third subset and hence the 
overall redundancy will be reduced by one. The idea is to use 
the technique of subsection lVI-AI The parity-check matrix Hi, 
in subsection IVI-AI was used to find the pattern of the error, 
and the parity-check matrix H 2 was used to find the location 
of a burst given its pattern up to a cyclic shift. 

This technique can be simply generalized for two- 
dimensional and multidimensional codes. We will describe it 
only for two-dimensional codes. Each of our constructions for 
two-dimensional codes uses two components codes Ci and 
C 2 . Assume that Ci and C 2 are binary codes with n x m 
parity-check matrix Hi and r 2 x n 2 parity-check matrix 
H 2 , respectively. We construct a two-dimensional parity-check 
matrix H for our two-dimensional code as follows. H is 
a two-dimensional matrix whose shape is the shape of the 
two-dimensional codeword. In each position T in this two- 
dimensional shape, TL has a column vector of length r\ + r 2 
which is a concatenation of a column from Hi and a column 
from H 2 . The column from Hi is £i if £i is the color given 
to position T by the first coloring (as mentioned before all 
our constructions can be represented by the coloring method). 
Similarly, the column from H 2 is £2 if £2 is the color given to 
position T by the second coloring (to avoid confusion, in the 
coloring method all positions are assigned with a color, by 
each coloring, including the redundancy bits). Now, we can 
use a proof similar to the proof of Lemma [13] to show that TL 
is a parity-check matrix for the required code. A generalization 
for multidimensional codes is straightforward. 

Finally, note that the same method can be also applied to 
the three constructions presented in [8]. Hence, we can supply 
a parity-check matrix for each code constructed by these three 
constructions. 

IX. Conclusion and Open Problems 

As we wrote in the abstract, the main results of the paper 
are summarized as follows: 

1) A construction of small redundancy multidimensional 
codes capable to correct a box-error. These codes and the 
box-error have considerably more flexible parameters 
from previously known constructions. 

2) A novel method based on D colorings of the D- 
dimensional space for constructing Z?-dimensional codes 
correcting a D-dimensional cluster-error of various 
shapes. 
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3) A transformation of the £>-dimensional space into 
another D-dimensional space in a way that a D- 
dimensional Lee sphere is transformed into a shape 
located in a £>-dimensional box of a relatively small 
size. This transformation enables us to use the previous 
constructions to correct a Z?-dimensional error whose 
shape is a D-dimensional Lee sphere. 

4) Applying the coloring method to correct more efficiently 
a two-dimensional error whose shape is a Lee sphere. 

5) A construction of one-dimensional and multidimensional 
codes capable to correct a burst-error of length b in 
which the number of erroneous positions is t. 

6) Applying the construction for correction of a Lee sphere 
error and the construction for correction of a cluster- 
error with small number of erroneous positions, to 
correct a Z?-dimensional arbitrary cluster-error. 

All the codes we have constructed are binary. We didn't 
discuss cluster-correcting codes over GF(q), but most of our 
results can be generalized straightforward for codes over 
GF(q). We have omitted some tedious proofs. The interested 
reader is referred to [18] to see some of these proofs. 

Clearly, our constructions do not cover all possible pa- 
rameters. Moreover, the redundancy of our codes is close to 
optimal, but not optimal, so there is lot of ground for possible 
improvements with possibly new construction methods. In 
fact, the main disadvantage of our methods is that, for large b, 
the lengths of the known &-burst-correcting codes and &-burst- 
locator codes are very large. 

Another question we didn't discuss in this paper is con- 
structions for cluster-correcting codes which correct a small 
cluster. This question was considered in [17]. The construction 
of D-dimensional cluster-correcting code which corrects a D- 
dimensional Lee sphere error with radius one is important in 
this connection. As we mentioned in subsection IV-CI we don't 
know the maximum number of (2D + l)-burst-locator codes, 
among the D component codes, that we can use. 

The next question is how to implement the coloring method 
for correction of multidimensional Lee sphere errors with 
dimension greater than two and radius greater than one? 

Finally, we still don't know and even don't have any indica- 
tion what should be the excess redundancy of an optimum code 
which corrects an arbitrary multidimensional cluster-error. 
Even the two-dimensional case is far from being resolved. 
The gap between the lower and upper bounds of Theorem [9] 
and Theorem [8] respectively, is quite large and we believe that 
both bounds can be improved. 
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